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^3 ■ 

■ We derive inversion formulas involving orthogonal polynomials which can be used 
^ ' to find coefficients of differential equations satisfied by certain generalizations of 

' the classical orthogonal polynomials. As an example we consider special symmetric 

^■f^ ^ generalizations of the Hermite polynomials. 

' 1 Introduction 

■ In Q we found differential equations of spectral type satisfied by the generalized 
\ Laguerre polynomials which are orthogonal on the interval 

■ [0,oo) with respect to the weight function 

-a;"e-^ + M6{x), a > -1, M > 0. 

T{a + 1) ^ ^ 

. , These orthogonal polynomials were introduced by T.H. Koornwinder in 0. 

' In order to find the coefhcients of these differential equations we had to solve 

, systems of equations of the form 

T~H ' oo 

g: ^a,(:r)i?%(r)(x)=K(x), n= 1,2,3,..., (1) 

d U 

5^ , where D — — denotes the differentiation operator. In tJ H. Bavinck showed 

dx ^ 

that the coefficients {ai{x)}'^^ are uniquely determined and can be written 



' in the form 

s 

X 



(2;)^(-ir^i|~^)(-2:)F,(x), z = l,2,3.... (2) 



This result is based on the inversion formula 



E ^tr'"'^(-^)4-7V) ^S,,, J <i,i, J = 0,1,2,.... (3) 



See also □. This inversion formula was derived in a similar way as t 
version formula involving Charlier polynomials found in a. See also 



,e m- 
and 
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section 3 of this paper. For inversion £Drinu]|as involving Meixner polynomials 
the reader is referred to a. See also □. In til we used the inversion formula 
to find differential equations of spectral type satisfied by the Sobolev-type 
Laguerre polynomials which are othogonal with respect to 

the Sobolev-type inner product 



T^T^TT / x"e-y{x)g{x)dx + M/(0)g(0) + Nf'{Q)g'{Q), 
r(a + 1) Jo 

where qlj> -1, M > and iV > 0. 

In Ej T.H. Koornwinder also introduced the generalized Jacobi polyno- 
mials {Pn'^''^^'^ (x)} which are orthogonal on the interval [—1,1] with 
respect to the weight function 

where a > —1, (3 > —1, M > and iV > 0. In we were looking for 
differential equations of spectral type satisfied by these generalized JLacobi 
polynomials. The general case turned out to be very difficult, but in 113 we 
were able to solve this problem in the special case that /3 — a and N — M. 
In order to find the coefficients of these differential equations we had to solve 
systems of equations of the form 

oc 

J2c^{x)D'P!{'■'^Hx)^Pn{x), n = 1,2,3,.... 

In i we showed that the coefficients {ci{x)}°^i are unique and that they can 
be written in the form 



This result is based on the inversion formula 
^ a + /3 + 2fc + 1 

k=i 

xPtr"'^-'—''ix)Pt^'''^'^(x)^S.,, J<^, z,j-0,l,2,..., (4) 

which is proved in i. This inversion formula was derived in a completely 
different way than the inversion formulas mentioned before. In u it is shown 
that this injzersion formula (with P — a) |-can be used to derive the results 
obtained in Ej in an easier way. Finally in E2I this inversion formula is used to 
solve the problem for all a > — 1, /? > — 1, A/ > and > 0. 



a + (3 + 2j + 1 {-a-i-i.-/3-i-i). 



^ {a + p + k + j + 
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In this paper we will derive several kinds of inversion formulas and we 
will show how they can be applied to find coefRcients of differential equations 
for generalizations of some classical orthogonal polynomials. 



2 Some classical orthogonal polynomials 

In this section we will recall some formulas involving classical orthogonal 
pobmomials which wc will use in this paper. For details the reader is referred 

toM. 

The Mcixncr-PoUaczek polynomials \ Pn'^\x; can be defined by 

their generating function 



(1 _ e^'At)-'+" (1 - e-^t)~'"" = £ pW(x; 0)^. (5) 

n=Q 



The classical Jacobi polynomials I P,["'^\x) \ can be defined for all a 
and P and n S {0, 1,2,.. .} by 

p(a,/3)^ ^ _ ^ (n + a + /3 + l)fc (a + fc + l)„^fc f x ~ 1 ^ 
They satisfy the orthogonality relation 

2"+'3+ir(a + l)r(/3 + l)y_i^ ^Mi + 2;j ^™ lxJ/^„ (x)dx 
a + + 1 (- + lW/3+l)n.,_^,„ = o,l,2,.... 



2n + a + /3+l (a + /3 + l)„n! 

The Gegenbauer or ultraspherical polynomials I Gn'^' (a;) > form a special 
case of the classical Jacobi polynomials. In fact we have 

Gi'H^) - -^^Pf A > -i, A ^ 0. (7) 

(A + 2)n ^ 

These ultraspherical polynomials can also be defined by their generating func- 
tion 

oc 

(l-2xi + t2)-A^^G(^)(^)r. (8) 

n=0 
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The special case A = needs another norniahzation. In that case we have the 
Chebyshev polynomials of the first kind {Tn{x)}1^_Q given by 



Tnix) 



P 



Pn " ''{1) 



2F1 



-n, n 
1 
2 



1-x 



n = 0,1,2,.. 



Their generating function equals 

1-xt 



l-2xt + i2 



J2Tn{x)e. 



(9) 



n=0 



The Chebyshev polynomials of the second kind {Un{x)}'^^Q are given by 

1-x^ 



Unix) = (n+1) " ' = (n+1) 2F, 



-n,n + 2 
3 
2 



, n = 0,l,2,.... 



These polynomials can also be defined by their generating function 

1 



l-2xt + i2 



(10) 



n=0 



Finally the classical Legendre (or spherical) polynomials {Pn{x)}^^Q form 
another special case of the classical Jacobi polynomials. In fact we have 



p„(x)=Fr'(^)=E 



(n + fc)! f X — 1 



k=0 



(n- fc)!(fc!)2 



n = 0,l,2,. 



These Legendre polynomials can also be defined by their generating function 

1 



Vl - 2xt + i2 



(11) 



n=0 



Note that the Legendre polynomials also form a special case of the ultras- 
pherical polynomials, since we have 

Pnix)=Gi^\x), n = 0,l,2,.... (12) 

The classical Laguerre polynomials can be defined for all 

I- J ra=0 

a and n G {0, 1,2,.. .} as 



4"w^D-i)'(:!:)i;^E(-i) 

fe=o ^ ^ fe=0 



. {a + L: + l)a^k x^ 
(n-fc)! fc!' 
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The generating function for the classical Laguerre polynomials is given by 

xt 



lexp — -U5]4")(x)r. (13) 

Further we have for n = 0, 1, 2, . . . 

i^^4")(a;) = (-ir4"_V)(x), z = 0,l,2,...,n. (14) 

Another family of continuous orthogonal polynomials is the one named 
after Hermite. The classical Hermite polynomials {Hn{x)}'^^Q can be defined 
by their generating function 



exp fxt - ^tA = ^ Hn{x)e. (15) 

Here we used another normalization than in 0. This one turns out to be 
more convenient in this paper. These classical Hermite polynomials satisfy 
the orthogonality relation 

— / e''' H„,{x)Hn{x)dx ^ m,n = 0,1,2, ... . 

J-oo ^ 

Further we have for n = 0, 1, 2, . . . 

D'H^{x)=H^-^{x), i = 0,1,2,..., n. (16) 
From the generating function it follows that 

if2„+i(0) = and i72„(0) = n = 0,l,2,.... (17) 

Further wc will use the kernels 

if„(a;,y) =^2*'-fc!i7fc(a;)iJfe(2;), n = 0,l,2,.... (18) 

fe=0 

By using ( p7| ) we easily find that for n = 0, 1, 2, . . . 

if2n+i(x, 0) = K^^{x, 0) = Y.^~\f^^H2k{x) (19) 



and 



K^^UO, 0) = i^.„(0, 0)^± = E % = (20) 

fc=0 ^ '' k=0 
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Finally we will consider the discrete orthogonal polynomials named after 
Meixner and Charlier. We choose normalizations different from those in tJ. 
The classical Meixner polynomials ■! Mn^\x] c) \ can be defined by their 

L J n=0 

generating function 

■.X CO 

1-- (l-<)-^~^ = EA4«(x;c)r. (21) 

The Meixner polynomials are connected to the classical Jacobi polynomials 
in the following way 



M^f\x; c) = P^/3-i.--/5--) , n = 0, 1, 2, , 



(22) 



The classical Charlier polynomials | cli''^ (x) \ can also be defined by their 

L J n=0 

generating function 

oo 

e-^^l+tf ^Y.Ci^\x)r. (23) 



3 Some inversion formulas 

In i we observed that the generating function (|2^) implies that 

oo / n \ 
n=0 \fc=0 / 

which implies that 



A:=0 

or 



I — n ^ 



1, n = 

0, n= 1,2,3... 



E {x) = S,, , j < i, z, j = 0, 1, 2, . . . . (24) 

k=j 



As already indicated in □ this formula (24) can be interpreted as follows. If 
we define the matrix T — {tij)2j=o with entries 

10, 3> i, 
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then this matrix T is a triangular matrix with determinant 1 and the inverse 
f/ of T is given by T^^ = U = {uij)2j=Q with entries 

_/ct;)(-.T),j<z 



0, j > i. 

Therefore we cah ( p^ ) an inversion formula. 

In the same way we find from the generating function (|l^) for the classical 
Laguerre polynomials 

i 

J2 LtV\-x)L[%{x) = J < J = 0, 1, 2, . . . . 

However, in view of ( p^ ) this inversion formula cannot be used to solve the 
systems of the equations of the form (|l|). In El H. Bavinck observed that it 
also follows from the generating function ( [T^ ) that 

(1 - ty-'-' = (1 - t)---'~' exp ( (1 - tr+^ exp ^ -""^ 



J-lJ \t-l 

oo / n \ 

which implies, by comparing the coefHcients of <*^^ on both sides, that 
g Lt^''^(x)Lt^:r'\-x) = J < z, z, J = 0, 1, 2, . . . , 

k=Q 

which is equivalent to This inversion formula implies that the system of 
equations ([l]) has the unique solution given by (|^). 

4 More (inversion) formulas 

Applying the method described in the preceding section to the generating 
function ( |lo|) for the Chebyshev polynomials of the second kind and the gen- 
erating function ( pT|) for the Legendre polynomials we obtain 

1 1 



oo 



n=0 



OO OO oo / n \ 

^Pfc(x)t^- ^ Pm{x)t"' = E [Y.Pk{x)Pn-k{x) f 
k=0 m=0 n=0 \k=0 ) 
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which imphes that 

n 

Y,Pk{^)Pn-k{x)^U„{x), n = 0,1,2,. 



k=0 



Another interesting formula of this kind can be found by using the generating 
function (|lO|) for the Chebyshev polynomials of the second kind and the gen- 
erating function for the Chebyshev polynomials of the first kind. In fact, 
we have 



oo 



1 1 - 



71=0 



2xt + t2 1-xt VI - 2xt + t2 



fc=0 m=0 n=0 \k=0 / 

which implies that 



J2 x'^T^-kix) = C/„(x), 71 = 0, 1, 2, ... . (25) 

k=0 

As before we can use the generating function (|^) for the ultraspherical 
polynomials to obtain the inversion formula 

■i 

^ Gtl\x)G^^l^(x) = J < z, z, J = 0, 1, 2, . . . . (26) 

k=j 

In view of ( p^ the special (limit) case A = ^ should lead to an inversion 
formula for the Legendre polynomials. If we define for every positive integer 
TV the matrix A — (ay )f^-^]^ with entries 

_\P,.,{x),3<i 

0, i > 



then this matrix is a triangular matrix with determinant 1 and hence invert 
ible. Now we have G^^^ {£) — 1, G^^\x) — 2Xx —x for A ^ — ^ and for 



n = 2,3,4,... 

Gi'Hx) = -^^Pi'-^-'-^\x) ^ ^Pi-'^-'\x) for A^-i. 

(A + i)„ n-1 2 

Now we have by using (^) for n = 2, 3, 4, . . . 

S„(x) := _^p(-i-i)(a:) = 1(1 - x)Pi'_:-'\x). (27) 
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bij 



Hence, the inverse A ^ = B = is given by 

f 0, i<j 

1, i = j 

—X i = j + 1 

In case of the Chebyshev polynomials of the second kind we can obtain 
an inversion formula as follows. If we define for every positive integer A'" the 
matrix A = {aij)fj^i with entries 



_ jUi-j{x), j < i 



then this matrix is a triangular matrix with determinant 1 and hence invert- 
ible. It is not difficult to show that its inverse A~-^ = B = {bij)fj^i is given 

by 

1, i = j 
—2x, i = j + 1 
1, i=j + 2 
0, otherwise. 



This can be shown by writing 



N 



BA = C = {cij)^j^-^ with Cij = '^bikakj 



k=l 



and showing that C 
well-known relation 



/, the identity matrix. This is done by using the 



Unix) - 2xUn+i{x) + Un+2{X) = 0, n = 0, 1, 2, . . . . 

In case of the Chebyshev polynomials of the first kind we consider the 
matrix A = {aij)fj_i for every positive integer N with entries 



■jix),j<i 
j > i. 



Then this matrix is also a triangular matrix with determinant 1 and hence 
invertible. The inverse A~^ = B = is given by 



bu 



fo, 

1, 



-X, 



I < J 
i=j 
i = j + l 



v'-^-'^{l-x'^),i>j + 2. 
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This can also be shown by writing 

N 

BA C = {cij)^j=i with = bikUkj 

fc=i 

and showing that C = I, the identity matrix. This is done by using the 
formula (25) and the well-known relation 

{l~x^)Un{x)-xTn+i{x)+Tn+2{x)=0, n = 0,1,2,.... 

The generating function can also be used to find inversion formulas 
involving Meixner-PoUaczek polynomials. In fact we have 

i 

J2 (-^; (2^; 0) = <5,, , j < j = 0, 1, 2, . . . 

k=i 

or 

i 

J2 Ptk ~<f)Pl^^, </>) = % , J < J = 0, 1, 2, . . . . 

By using the generating function ( pT|) for the Meixner polynomials we 
find the inversion formula 

i 

J2 Mtk^ i-x; c)Mi^}^ {x- c) = % , J < z, J = 0, 1, 2, . . . . (28) 

k=3 

We remark that this inversion formula is different from the one obtained ini. 
See also □ for an application of that inversion formula. 

Note that the generating function ( [l5| ) for the classical Hermite polyno- 
mials implies that 

1 = exp (xt~-A exp (-xt+h^\ = ( i?fc(a;)-ff„-fc(ix)z"-M 
which implies that 



^fffe(x)iJ„_fc(ia;)r-^- - <^ 
fc=o 



1, n = 

0, n 1,2,3. 



In view of ( |16D this formula can be used as follows. A system of equations of 
the form 

oo 

F^x) = akix)D''Hnix), n = 1, 2, 3, . . . , (29) 

fc=i 
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where the coefficients {«fc(a;)}^2 ^^'^ polynomials which are independent of 
n, has the unique solution 

k 

afe(a;) = ^i*'--'fffe_,(ia;)^;(a;), fc= 1,2,3,.... (30) 

5 Inversion formulas involving Jacobi polynomials 

In i we have found the inversion formula involving Jacobi polynomials. 
As mentioned before this formula was found in a completely different way. 
The well-known generating function for the classical Jacobi polynomials has 
a different structure, pwhich implies that the method used before cannot be 
used in that case. In u we proved that for n = 0, 1, 2, ... we have 

Now y — X leads to the inversion formula (^. If y — —x this leads to a 
formula which was used in □ in the case that (3 = a. 

By using the relation (|2^) between the Meixner and the Jacobi polyno- 
mials we find from the inversion formula (28) for the Meixner polynomials 
that 

E Ptr"'''''^'\^)Ptt'^'\^) = J < z, J = 0, 1, 2, . . . . 

Another inversion formula involving Jacobi polynomials can be obtained 
from the inversion formula (^6|) for the ultraspherical polynomials. By using 
(0) and after setting \ — a +^ this leads to 

i 

E 



( -2a-l),_fc (2a+l)fc_j 



X Ptk~''~''~'\^)Pi-fi^) - 3 < h hJ = 0, 1,2, . . .. 
6 Applications to diflferential equations 

In this section we will investigate the generalized Hermite polynomials 
{iJ^(a;)}^_g which are orthogonal on the real line with respect to the weight 
function 

w{x) = -^e"-"^ + MS{x), M > 0. 
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In □ these generalized Hermite polynomials are called special (linear) per- 
turbations of the classical Hermite polynomials. They can be represented in 
terms of the kernels (|l^) as (see □) 

H^{x) = Hn{x) + A/g„(x), n = 0, 1, 2, . . . , 

where Qq{x) = and 



H„{x) Kn-iix,0) 
H„{0) i^„-i(0,0) 



^qn.kHkix), n = 1,2,3, 



k=0 



with, by using <^^), for n = 1,2,3,... 

g„,„ = i^„_i(0,0) and = -2'= fc! i/fe(0)H„(0), fc - 0, 1, 2, . . . , n - 1. 

In i it is shown that these generalized Hermite polynomials satisfy a differen- 
tial equation of the form 

oo 

M ak{x)y^''\x) + y"{x) - 2xy'{x) + {2n + Ma„)y(a;) = 0, (31) 
fe=i 

where the coefficients {afc(x)}^j^ are polynomials with degree [afc(x)] < k, k = 
1,2,3,... which are independent of n. Moreover it is shown that the 'eigen- 
value' parameters {a2ra+i}^o '^^^ ^® chosen arbitrarily, 

n 

ao = and a2n = '^{^2] - ^2^-2)92^,2^, n = 1, 2, 3 . . . , 

i=i 

where A„ = 2n, n = 0, 1, 2, . . .. Hence, X2j — )^2j-2 =4, j = 1, 2, 3, . . . and by 
using (@) 



q2j,2j = i^2,-i(0, 0) = i^^, J - 1, 2, 3, ... . 



-1 

(j - 1) ' 
This implies that 



" (3\ (3\ /'5^ 



^(7-1)! ^ fc! (n-1)!' 

In order to find the coefhcients {Qfc(x)}^i we set j/(a;) = H^{x) = Hn{x) + 
MQn{x) in the differential equation ( |3l| ) and view the left-hand side as a 
polynomial in M . Then the coefficients of this polynomial must vanish, hence 

00 

J2ak{x)D''H„{x) ^ ~ar,H„{x)-Q';{x) + 2xQ',^{x)-2nQ„{x), n = 0,l,2,... 
fc=i 
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and 

CXD 

^ ak{x)D''Qn{x) = -a„Q„(x), n = 0, 1, 2, . . 
fc=i 

Since we have, by using and (po|), 



(S.) 

Q2n+l{x) ^ K2n{0,0)H2n+l{x) = ^H2n+l{x), n = 0,1,2, 



l2_ 

n< 

and 



%^0, n = 0,l,2,... 
both systems of equations lead to 

2ji+l 

afe(x)if2ri+i-fc(a;) = -a2n+i^2«+i(a;), n = 0,1,2,.... (32) 

k=l 

Further we have 

n 

Q2nix) = ^g2«,2feff2fe(a;), n = 1,2, 3, . . . , 
where, by using (|2C|), 

g2n,2n = -fSr2n-l (0, 0) = , ^ "T,, , n = 1,2,3, .. . 

(n — Ij! 

and, by using (p7|), for fc = 0, 1, 2, . . . , rt — 1 and n = 1, 2, 3, . . . 

g2„,2fc = -22^-(2fc)!ff2fc(0)i/2„(0) = ^"^^irjlf^^' - 
Now we have for n ^ 1,2,3,... 



Q'i,ix) - 2xQ'2„{x) = ^g2„,2fe [H!;,{x) - 2xff^,(a 

n 

= ^ 92n.2fe [-4kH2k{x)] . 

Hence, for n = 1, 2, 3, . . . we obtain 

-Ol2nH2n{x) - Q'^ni^) + 2xQ'2„ix) ~ 'inQ2n{x) 

n 

= -~a2nH2n{x) - 4^(n - k)q2nakH2k{x), 

k=0 
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which leads for n = 1, 2, 3, ... to 



2n n 

'^ak{x)H2n-k{x) = -a2nH2n{x) -4^(n- k)q2n.2kH2k{x). (33) 
k=l k=0 

Hence, with and (|3^) we have found that 

oo n 

J2ak{^)D''Hn{x) ^J2''k{x)Hn-k{x) ^ Fn{x), 1,2,3,..., 

fe=l fe=l 

where 

-F2„+i(a;) = -a2„+ii/2n+i(a;), n = 0, 1,2, . . . 

n 

F2nix) = -a2nH2n{x) - 4^(n - k)q2„,2kH2k{x) , n = 1,2,3, ... . 

fc=0 

This system of equations has the form (|29|). So we can use ( |30| ) to conclude 
that 

fe 

afc(a;) = ^i^-^i/fc_,(ix)^;(a;), fc= 1,2,3,.... 

We emphasize that these generahzed Hermite polynomials are orthogonal 
with respect to a weight function consisting of the classical Hermite weight 
function and a Dirac delta distribution at the origin. Therefore these general- 
ized Hermite polynomials could be considered as Krall-Hermite polynomials, 
but, these are quite different from the Krall-Hermite polynomials considered 
in □ which are pot orthogonal. 

Finally, in t3 it is shown that these generalized Hermite polynomials can- 
not satisfy a finite order differential equation of the form (pl|). 
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